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DEVELOPMENTS IN THE SIMULATION OF 
COMPRESSIBLE INVISCID AND VISCOUS FLOW 
ON SUPERCOMPUTERS 

Joseph L. Steger and Pieter G. Buning 
INTRODUCTION 

The near term availability of scientific supercomputers will soon permit routine 
simulation of three dimensional compressible flow about relatively complex config- 
urations. 

General purpose flow simulation codes have been under development at NASA 
Ames for application to existing large scale computers (such as the CRAY-XMP 
and the Cyber 205) and future supercomputers. These codes have primarily used 
the finite difference or finite volume approach, but because of storage and speed 
limitations, the codes have previously been restricted to simulation of fairly simple 
geometric configurations. However, work towards generalizing the current codes to 
treat complex geometries and larger sized meshes has been underway for the past 
several years in anticipation of supercomputers. 

While there is no agreement within NASA Ames as to what is the best general 
purpose simulation procedure, considerable experience has been obtained with a 
class of implicit finite difference algorithms written is terms of generalized coordi- 
nates. These implicit procedures are widely used in various CFD simulation codes, 
and they are rapidly being extended to treat quite complex geometries using cur- 
rent and future large scale computers. This paper will review these procedures 
and discuss how they are being deployed on supercomputers for simulation of com- 
plex three dimensional flow fields by using various mesh interface schemes. The 
importance of flow visualization and diagnostic methods to three dimensional flow 
simulation will be discussed at the end of this paper. 

FINITE DIFFERENCE ALGORITHMS 

A general purpose simulation code should be able to solve either steady or un- 
steady, viscous or inviscid flow with only minor input changes from the user. While 
it is unlikely to be an optimum code for all tasks, it should be good at its mainstay 
tasks, and have acceptable efficiency for other tasks. Ideally it will be a modular 
code whose basic engine can be stripped out and readily applied to other problems. 
Our mainstay general purpose flow simulation code has been a centrally differenced 
implicit approximate factorization code [1-6] which solves a conservative form of ei- 
ther the Euler or thin layer Navier Stokes equations cast in generalized coordinates. 



1. Fluid Conservation Equations and Transformations 

The conservation equations of mass momentum and energy referenced to a Carte- 
sian coordinate system can be represented in the flux vector form (c.f. [7]): 

d t Q + d x {F+ F v ) + d y {G + G„) + <M H + H V )= 0 (1) 

The viscous flux terms F v , G v and H v contain derivatives and throughout a nondi- 
mensional form of the equations will be used. The conservative form of the equa- 
tions is maintained chiefly to capture the Rankine Hugoniot shock jump relations 
as accurately as possible. 

New independent variables r, £, »/, f are generally chosen to map a curvilinear 
body conforming discretization into a uniform computational space as shown in Fig. 
1. Body conforming curvilinear meshes are generally used in finite difference and 
finite volume computations for a variety of reasons: to simplify the application of 
boundary conditions, to allow clustering of grid points to flow field action regions, 
to help maintain the well-orderliness that is useful for vector processing and for 
various implicit methods that employ approximate factorization techniques such as 
ADI, etc. Indeed, the calculation of high Reynolds number viscous flow is simply 
impractical without the use of a body conforming curvilinear mesh that employs 
clustering in the direction normal to the body surface. 

In the new independent variables the transformed equations can be represented 
as (c.f. (7] for the detailed terms): 

d T Q + d ( {F + F v )+d 1 ,(G + G v ) + df[H + H v )= 0 (2) 

where the original dependent variables are maintained. (The flux vectors of the 
transformed equations can be made to resemble their Cartesian counterparts by 
combining terms into contravarient velocity components [3-4].) If a body conforming 
coordinate is used, then for high Reynolds number flow it is generally permissable 
to make a thin layer assumption and to discard viscous terms except for those in 
the normal-like direction. If f is the coordinate away from the surface, the thin 
layer equations can be represented as [3,4] 

d T Q+d ( F + d n G + d ( H =Re~ l d<S (3) 

where the viscous terms in f have been collected into the vector S and the nondi- 
mensional reciprocal Reynolds number is extracted to indicate a viscous flux term. 

In differencing these equations it is often advantageous to difference about a base 
solution denoted by subscript o as 


Sr{Q ~ Qo) +8t(E- E 0 ) + S V (F - F 0 ) + 8 ( {G- G 0 ) - Re~ l 8 ( (5 - S 0 ) 

= —d r Qo — d{Eo — dqFo — d t Go + Re l d ( So 

where 8 indicates a general difference operator, and d is the differential operator. 
If the base state is properly chosen, the differenced quantities can have smaller 
and smoother variation and therefore less differencing error. If a the base solution 
exactly satisfies the partial differential equation, then the right hand side of Eq.(4) 
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is identically zero. A uniform free stream satisfies the governing equations and 
is frequently taken as the base solution so as to minimize far field differencing 
errors. Such an error can occur if the coordinate transformation terms that have 
been embedded into the transformed fluxes are not consistently differenced [3,4,8]- 
There are also anologies between the finite volume method and the transformed 
difference equations, and the actual discretized equations that result from the two 
schemes can be made similar. 

2. Implicit Central Difference Algorithm 

Equations (3) or (4) have been solved using a Beam- Warming noniterative ap- 
proximate factorization implicit scheme of the form [3-6] 

[j + W e A"-Di| e ] [l + h6 n B n -D,|„] + A Q" 

= —At [*,(£" -E 00 ) + 6 n (F n -F 00 )+6 i (G n -G 00 )-Re- l 6 s (S n - S TO )] (5) 

-D c (Q n -Qoo) 

where h = At or (At)/2 and the free stream base solution is used. Here S is 
typically a three point second order accurate central difference operator, while 6 is 

A A, 

a midpoint operator used with the viscous terms. The matrices A,B,C, and M 
result from local linearization about the previous time level and J is the Jacobian 
of the coordinate transformation. The factored left hand side operators form block 
tridiagonal matrices. 

Because central space difference operators are used, numerical dissipation terms 
denoted as Di and D e have been inserted into Eq.(5). In their simplest form these 
have been given as combinations of fourth differences 

D e = e e At J-M(VA) 2 e + (VA)J + (VA)J) J (6a) 

and second differences 


A|{=«(A(J’- 1 (VA) { J, D,- 1, =s,A(/-'(Vd),J ! A| t =£,a;j-‘(va),j 

(6b) 

where V and A are twopoint backward and forward difference operators and where 
e,- > 2e tf , and e e = O(l). The implicit second difference numerical dissipation 
operators were chosen to keep the left hand side factors block tridiagonal, and as 
Di works on A Q, accuracy is not impaired. More robust dissipation terms are 
often used in which the e coefficients are scaled with some approximate modulus of 
the A, B, and C Jacobian matrices and gradients of an appropriate variable. For 
example, explicit and implicit ^-dissipation terms such as 

D e = (Atf)J -1 (|£ r | + |£ y | + l£*D[ e2 ^|(J + + 64< ^ J 

Di = (Ai) J -1 (|fj| + |£ y | + 8 + 3 e 4 P]J 

have been used where p is the nondimensionalized fluid pressure and e 2 is 0(1) while 
€4 is 0(0.1). In transonic flow the term (|^| + |$ y | + |£ z |) is an estimate to the 
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spectral radius of A for variables that have been nondimensionalized with respect 
to the sound speed. Since (|£ 3 | + |£j,| + |£ z |) does not contain a fluid variable, it 
can be left as an outside coefficient to the dissipation operator without effecting the 
weak conservation form of the difference equations. 

Body surface boundary conditions have usually been supplied by using a combi- 
nation of normal-momentum, tangency or no slip, and extrapolation [3,4]. Various 
far field conditions have been used including characteristic-like conditions [9-11]. 
Because of their simplicity, in most of the application codes the boundary condi- 
tions have been imposed explicitly, or a combination of simplified implicit-explicit 
conditions [12] have been used. However, fully implicit boundary conditions have 
also been used ( c.f. [1,11,13]) and an elegant implicit characteristic-like formulation 
has been given and tested by Chakravarthy[ll]. 

3. Vectorization and Multi-Tasking 

The structure of the above algorithm lends itself to vectorized computer coding. 
Vectorization has been implemented by inverting ’’pencils” of data [14], but FOR- 
TRAN based codes generally follow an approach coded by Benek (unpublished, 
circa 1980). In this approach strings of block tridiagonal matrices are inverted si- 
multaneously so as to avoid the recursive nature of matrix elimination procedures. 
For example, if indices j, k, l correspond to £, q, f, the first factor of Eq.(5) forms a 
block tridiagonal between, say, points j = 1 to j = jmaz. There is one such £-block 
tridiagonal for each k, l index so by simultaneously inverting £-block tridiagonals 
over, say, k = 1, kmax a vector length of kmax is achieved, see Fig. 2. However, 
in inverting a £-block tridiagonal temporary storage is needed for the backward 
elimination and this storage requirement is increased by kmax. (Means of reducing 
the block size are discussed in [15-17] as a way to improve efficiency, and these same 
techniques reduce temporary storage as well.) Overall, vectorization tends to com- 
plicate the coding, but not unduely so. On the CRAY-XMP the ARC2D vectorized 
code [5] runs about 5 times faster than the unvectorized code for a two-dimensional 
grid of order 200 x 40. 

The central differenced implicit algorithm, Eq.(5), has also been experimentally 
coded for multi-tasking on two minicomputers by assigning each processor a portion 
of the code [18]. For example, processor 1 can solve the £-block tridiagonal in the 
range l = l,lx, processor 2 can solve the £-block tridiagonal from l = li + 1, l 2 and 
so on. Because of a lack of readily available hardware and software, multi-tasking 
has not yet worked its way into the codes used for routine applications. 

4. Improvements in Efficiency and Accuracy 

The introduction of a faster computer (that is still easy to use) drives the develop- 
ment of much more efficient numerical algorithms. This is because as the machines 
become more powerful, it is easier to experiment will new ideas in numerical algo- 
rithms. New algorithms can be quickly verified, problem areas can be isolated, and 
optimum convergence parameters can be found. As a result, improvements in nu- 
merical algorithms have kept pace with improvements in computer hardware. Such 
has been the experience with the central differenced implicit algorithm. Originally 
coded for the Control Data 7600, it was transferred to the CRAY-1S and vectorized 
[5] to run an order of magnitude faster. Increased computer speed has in turn has- 
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tened the numerical optimization process indicated below which has lead to another 
order of magnitude improvement in the code’s steady state performance. 

Improvements to the basic algorithm in both efficiency and accuracy have been 
made by a variety of contributors. To improve its overall efficiency, simple changes 
have been optionally implemented into Eq.(5). These changes have been selected so 
that they do not unduely complicate the basic algorithm. They include the use of 
space varying At, use of a sequence of coarsened grids to provide a good initial guess, 
cutting inversion costs by using either diagonalization [15,16] or block reduction [17] 
methods, implementation of better numerical dissipation terms, and more implicit 
treatment of the numerical dissipation terms. As described in [5,6], these combined 
changes can improve steady state efficiency by an order of magnitude. Details and 
appropriate references are described in [5,6], 

Although the implicit algorithm has been presented with three point central dif- 
ferencing, versions of the algorithm that have fourth order accuracy in space have 
been available [1-4] and are preferred unless strong shock waves are present. Reddy 
[19] has also demonstrated a version in which a pseudo-spectral operator is used in 
place of the right hand side convection operators. Improved dissipation models [5,6] 
and total variation diminishing (TVD) implementations [20] have also been carried 
out to better capture shocks, and perturbation about approximate base solutions 
has been used to reduce the number of needed grid points [21]. 

5. Applications 

The factored implicit scheme described above has been ultilized in a variety of flow 
applications. Two-dimensional simulations include steady and unsteady flow about 
airfoils [c.f. 3,5,6,22-25], cascades [9], projectiles [26-27], and inlets[l6]. Three di- 
mensional calculations have been carried out about simple configurations using rel- 
atively coarse meshes. Simulations include the supersonic blunt body problem[28], 
simple bodies [4-6], wings [29], afterbodies [30], and the space shuttle [31]. Figures 
3 and 4 are reproduced from the publications of Diewert [30] and Chaussee et al 

[31] and show representative three dimensional viscous flow solutions that have been 
obtained on a CRAY-1S. To save computer time and storage, both of these calcu- 
lations where computed in stages. In the boattail case a complete solution was first 
obtained using a grid that inadequately resolved the base flow region. A refined grid 
was then introduced about the afterbody section and the computation for only this 
region was carried out using inflow conditions taken from the coarsened grid result. 
In the shuttle calculation a parabolized Navier Stokes code was used in those flow 
regions in which the main outer flow remains supersonic and the inner viscous layer 
remains attached with respect to the mainstream direction. These calculations and 
a majority of the viscous calculations carried out using the factored implicit algo- 
rithm at high Reynolds numbers have ultilized a simple algebraic turbulence model 

[32] , 

6. Upwind Schemes 

As an alternative to using central space differencing for the convection terms, 
upwind (i.e. backward and forward) space differencing can be used if the fluxes 
are properly split according to their characteristic properties. As discussed in [33] 
upwind schemes can have several advantages over central difference schemes, in- 
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eluding natural numerical dissipation, better explicit stability, and more readily 
inverted implicit schemes. Conversely, upwind schemes for systems of equations 
have generally been more complicated and computationally expensive than central 
difference schemes and are not very suitable for treating viscous terms. 

Although upwind schemes are not as extensively ultilized in aerodynamics simu- 
lations as central difference schemes, upwind schemes have been used on curvilinear 
grids in both finite difference and finite volume formulations and good results have 
been obtained. Figures 5a and 5b, for example, show an inviscid flow result for a 
NACA 0012 airfoil at M = 0.8 and an angle of attack of 1.25° using a flux split 
class of upwind scheme [34] (see [35] for similar results). Various upwind schemes 
have also proved to be quite effective for capturing strong shock waves, c.f. [36] 
and [37]. However, because of their complexity and cost, upwind based schemes 
have not yet been widely distributed as general purpose simulation codes for either 
inviscid or viscous flow. 

Upwind and partially upwind based schemes will likely first emerge as important 
flow simulation schemes in two areas: 1) supersonic and hypersonic flows in which 
there are very strong embedded shocks, and, 2) transonic and supersonic flows in 
which one coordinate direction can be mostly aligned with the flow streamlines. In 
this latter case, a mixed upwind-central difference scheme can be very effective if 
the upwind differencing is used in the streamline direction. For example, an implicit 
algorithm for the thin layer Navier Stokes equations could have the form 

[/+ h6 b ( {A + ) n + h6'C n - hRe~ 1 S { J~ 1 M n J - D,| f ] 

X [/ + A^(A-) n +A^5 n -Z7i| n ] AQ" = 

- A<{*‘[(£+)" - £+] + *'[(£-)" - E~] + 6 n (F n - £») 

+ 6 s (G n - Goo) - Re~ l S t (S n - Soo)} - D t {Q n - Qoo) (7) 

where and 5 ^ are backward and forward three-point difference operators and 
D e contains only // and f numerical dissipation operators. This two factor implicit 
scheme is readily vectorized or multi-tasked in planes of f = constant. A semi- 
implicit scheme is obtained by neglecting the calculation of h6^A~ in the implicit 

/N 

backsweep operating on A Q n . 

COMPOSITE GRIDS 

Use of a single well-ordered body conforming curvilinear mesh simplifies the ap- 
plication of boundary conditions and can lead to use of efficient solution procedures. 
However, the generation of body conforming well-clustered curvilinear grids that are 
not overly skewed and have smooth variation can often be quite difficult. In partic- 
ular, it is generally impractical to build a single grid of this type for complex three 
dimensional configurations. Of course, by judiciously introducing cuts in the grid, 
some fairly complex bodies can be meshed with a single grid, Fig. 6 illustrates this 
possiblility. But the trend has been to introduce more than one grid and to patch 
or overset the grid systems together. The sketches shown in Fig. 7 illustrate simple 
patch and overset grid configurations in two dimensions. (There has also been some 
activity in reverting back to the use of pure Cartesian grids that intersect the body 


6 



in a random way and so require special logic at the body interface [38]. Such a 
procedure appears to be impractical, however, for high Reynolds number viscous 
flow problems unless the Cartesian scheme is patched-to or overset-with a body 
conforming grid that is used only near the body. In this case the Cartesian mesh 
approach would be a special case of either the patched or overset grid method.) 

The use of a set of patched or overset grids to form a larger composite grid carries 
the discretization process one step further. In a sense such a finite difference pro- 
cess assumes some of the characteristics of a finite element scheme that uses large 
powerful elements in which each element is itself discretized. In this discretization 
process each individual grid in the system is well ordered and is thus suitable for ef- 
ficient finite difference solution using any available single grid scheme. The problem 
with such a composite grid scheme is the difficulty of interfacing each mesh without 
degrading numerical accuracy or convergence. Moreover, as few such meshes should 
be generated as necessary to achieve grid efficiency. 

Limited experience with both patched [39-41] and overset grids [42-46] has not 
shown which method is preferable - an optimum method will perhaps combine 
both patched and overset grids. Both schemes necessitate extensive bookkeeping 
procedures. The patched grid method has as its chief drawback a grid generation 
problem that is still relatively difficult because various interfaces have to be defined 
and grids have to be generated with both inner and outer defined boundary surfaces. 
Drawbacks to using overset grids include having to interpolate data points along 
an irregular boundary and the bookkeeping can be especially complex if more than 
two levels of overset grids intersect each other. 

Figures 8-9 show recently published two dimensional inviscid flow results ob- 
tained using patched and overset meshes. Both grids and flow contour levels are 
shown. The supersonic biplane results of Hessenius and Rai [41] were computed us- 
ing a patched grid method in which a flux balance interface scheme is used at grid 
boundaries to ensure fluid conservation. In this method [40], meshes must share a 
common boundary, but grid lines need not have a common slope or even join to- 
gether. These results were obtained using a first order accurate Osher scheme. The 
overset grid results of Dougherty [46] for an airfoil and flap were computed using 
the factored implicit algorithm. Only simple interpolation is used to interface the 
mesh boundaries using data from nearby points from the underlying grid. Figure 
9c shows the nearby interpolation points used to update the outer boundary of the 
flap grid (open symbols) and grid points which are excluded from the calculation 
(filled symbols). Simulation of wing-body configurations using the implicit algo- 
rithm with composite grids are currently underway by Holst and coworkers and 
Benek and coworkers and will appear shortly. 

POSTPROCESSING 


1. Flow Visualization 

One of the major problems that is being encountered in three dimensional flow 
simulation is the difficulty of displaying a limited and proper kind of data that 
will lead to better understanding. In the future, much of the computational aero- 
dynamicist’s time will be devoted to extracting and displaying various features of 
the solution. In three dimensions, flow phenomenon such as flow-reversals, shocks, 
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shear-layers, vortices, etc. can often be difficult to identify and visualize, especially 
if the flow is also unsteady. Graphic displays of contour surfaces and particle paths 
and the like can also be expensive to generate. As a result, more computer resource 
could be expended on analyzing a solution and displaying it in a meaningful way 
to a human being than what was needed to generate the solution in the first place. 
To become convinced of how costly this can become one need merely contemplate 
generating numerous detailed displays of a three dimensional flow field in which 
reliable hidden line removal is needed to keep the display from being too confusing. 

Two approaches to extracting graphic information and flow visualization are 
evolving. In one approach the engineer works at a graphics workstation and dis- 
plays various portions of the flow field, observing the solution from different vantage 
points, and otherwise interacts with the computed solution, perhaps by seeding par- 
ticles and observing their behaviour. Because the three dimensional data base can 
be quite large, this approach will require very powerful graphic work stations and 
high speed data links between the graphics station and the supercomputer. The 
advantage of this approach is that one can stumble across information that might 
not otherwise be anticipated. 

The other approach is to program the supercomputer itself to diagnose the data 
base. Such a computer program is needed in order to throughly search a flow field 
to find details that a human at a workstation would find too tedious to locate. 
Algorithms must be developed to search out and display special features such as 
shock waves, vortices, and separation lines. Because the data bases are so large, 
only the supercomputer itself will be able to perform many of these calculations. 

Reliable software to locate interesting flow regions is not available and will be dif- 
ficult to generate. For the CFD algorithm developer, developing flow visualization 
procedures offers interesting opportunities. Consider, for example, the difficulty 
of building an algorithm that automatically identifies embedded shock waves in 
a complex three dimensional flow, especially weak oblique shocks. Or consider 
the apparently straightforward task of building an efficient and accurate particle 
trace scheme. To trace particles given the velocity field, one need only numerically 
integrate simple ordinary differential equations. But if too many particles are dis- 
played, the picture will be confusing. Interpolation of the velocity components for 
each particle location can also be a costly process, so efficient algorithms are needed 
to minimize this expense. Finally, unless the particles are seeded in the right loca- 
tions, the most interesting features of the flow will not be observed. For example, 
particles seeded outside of a small vortex core will not be entrained inside the core. 

2. Diagnostics 

The accuracy of a finite difference solution is generally appraised by successively 
refining the grid and comparing solutions from one grid to another. If the solution 
is unchanged, the flow result is likely resolved. (Here we assume that a Reynolds av- 
eraged Navier Stokes solution can become invariant with grid refinement.) Because 
the next generation of supercomputers will not be ‘super enough’ to carry this pro- 
cess too far, the problem of determing solution accuracy may stimulate additional 
research in approximate methods such as potential and boundary layer schemes. 
We might compute, for example, a boundary layer solution to verify the correctness 
of the skin friction and heat transfer found from a Navier Stokes procedure. By 
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taking the Navier Stokes computed pressure gradient, the boundary layer equations 
can be solved using a very refined grid near the wall. If the boundary layer solu- 
tion returns the same wall values computed from the Navier Stokes equations, then 
the Navier Stokes solution is adequately resolved. (Because the boundary layer will 
likely be run in an inverse mode if the flow is separated, the check might be whether 
or not the boundary layer solution returns the same surface pressure distribution.) 

For complex three dimensional flows it will also be difficult to provide code veri- 
fication. Comparison with linear theory may be inadequate and experimental data 
may not be available or reliable. One possiblility is to compare the Navier Stokes 
computed results to those obtained from the solution of a simplified set of equations 
that are solved over the entire field, but in which rotational and perhaps even com- 
pressibility effects are directly taken from the Navier Stokes solution. As a quite 
simple illustration, the Poisson equation for the incompressible streamfunction tf> 

Ipxx + $ yy = -W 

should reproduce an incompressible two-dimensional Navier Stokes result if the 
vorticity, w, is taken directly from the Navier Stokes solution. Because this is 
a simple scalar equation, it should be possible to solve for it very quickly and 
using a more refined grid than what was used for the Navier Stokes simulation. 
If the refined grid Poissson calculation can return the Navier Stokes solution then 
considerable validation of the code is obtained. Carrying out such a process for 
complex configurations in three dimensions will be difficult but useful. 

CONCLUDING REMARKS 


Evolution of existing finite difference schemes should provide reliable and efficient 
single mesh codes for nonlinear three dimensional flow field simulation. To treat 
viscous flow about highly complex body configurations, however, composite meshes 
will be needed and considerable work in interfacing these schemes into multiple 
grid codes must yet be carried out. Just as the task of optimizing the flow solvers 
has been accelerated by the availability of large fast computers, so too the task of 
developing interface schemes will be accelerated with the availablity of forthcoming 
supercomputers, provided these machines are easy to use. 

Cost considerations are leading to supercomputers that achieve more and more 
of their performance advantage by use of repetitive processors. Many of these 
machines are more difficult to use than machines that achieve higher performance 
by having a faster single instruction time. Hopefully these future supercomputers 
will not be a machines of ‘insurmountable opportunity’. 

It is generally accepted within the CFD community that three dimensional data 
bases are so large that new ways of displaying and interacting with the data is nec- 
essary. It must be emphasized that hardware and graphics software provides only 
part of the solution. If the computed solutions are to be properly analyzed and 
understood, extensive new algorithm developments in llow visualization and diag- 
nostics will be needed, and these algorithms will require the same kind of creativity 
that went into generating the flow simulation procedures themselves. 
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Fig. 1. Sketch showing mapping of physical space to computational space for a 
well-ordered warped spherical grid. 
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Fig. 2. Simultaneous inversion of block tridiagonals to obtain a vector length. 



a). Streamline pattern for afterbody model 

Fig. 3. Transonic flow over a body of revolution with conical afterbody, at 8° angle 
of attack, = 0.9 and i?e<* = 3 x 10®. 
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Fig. 4. Various solution features of the space shuttle at = 7.9, 25° angle of 
attack and Rei = 6.07 x 10 s per inch . 

a) Surface streamline pattern. 

b) Cross sectional view showing grid slice conforming to outer bow shock. 

c) Pressure contours in vicinity of canopy region showing embedded bow shock. 

d) Cross sectional flow vectors. 


16 



- 1.5 



a) Surface pressure distribution, C p . 

Fig. 5. Calculation of inviscid transonic flow using flux vector splitting class of Suite 
difference scheme for NACA0012 airfoil at M a 0 = 0.8 and 1.25° angle of attack. 
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X 

b) Mach contours. 

Fig. 5. Concluded. 
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Fig. 6 Sketch illustrating the use of a single body-conforming warped-spherical grid 
with cuts. 
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Fig. 7. Sketches of a possible composite grid formed from patched and overset grids 
about a typical multiple body combination. 
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a) Mach contours 

Fig. 8. Inviscid supersonic flow about a biplane computed using a first order Osher 
scheme with patched grids, Moo — 1*5. 
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b) Composite grid formed using four patches. 


Fig. 8. Concluded. 
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a) Mach contours 

Fig. 9. Inviscid transonic flow about a generic airfoil detached flap combination 
using the implicit algorithm and overset grids, — 0.7. 
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b) Composite grid formed using two overset grids. 


Fig. 9. Continued. 
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c) Detail of minor grid showing all blanked points and interpolation points for 
the outer boundary of the minor grid. (Hole interpolation points not shown). 


Fig. 9. Concluded. 
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